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Abstract 



One relates factorization of bivariate polynomials to singularities of projective plane curves. 
One proves that adjoint polynomials permit to find the recombinations of factors modulo (x) 
induced by both absolute and rational factorizations, and so without using Hensel's lifting. One 
establishes in such a way the relations between the algorithms of Duval-Ragot based on locally 
constant functions and the algorithms of Cheze-Lecerf based on lifting and recombinations. One 
shows in such a way that a fast computation of adjoint polynomials leads to a fast factorization. 
The proof is based on cohomological sequences and residue theory. 
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1. Introduction 



Factorization of multivariate polynomials is a central topic in Computer Algebra. One 



refers to ( Cheze 2004 Cheze-Lecerf 2007 Gathen-Gerhard 2003 ) and to the references 



therein for recent surveys of the topic. In this article, one studies the relations between 
singularities of projective plane curves and factorization of bivariate polynomials by using 



adjoint polynomials. In (Duval 1991 ), an algorithm is given for absolute bivariate factor- 



ization based on locally constant rational functions on the curve, using normalization and 



rational Newton-Puiseux expansions. The best actual complexities for rational (Lecerf 



2007]) and absolute ( Cheze-Lecerf 2007 1 bivariate factorization have been obtained later 
on, based on a method of lifting and recombination of modular factors. One establishes 
here the bridge between these two approaches and one shows that the factorization can 
be computed fast from adjoint polynomials. 



Main result. Let F £ k[x,y] be a bivariate polynomial defined over a field k. We are 
interested in computing both the rational (over k) and absolute (over an algebraic closure 
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k) factorizations of F. In all of the sequel (except Section 7), one assumes that F satisfies 
the following hypothesis 

(H) F(0, y) is separable of degree d = deg(F). 

In particular F is square-free. Let C C P 2 be the (reduced) projective curve over k defined 
by F. One says that D C P 2 is an adjoint curve of C if it passes throw all singular points 
p of C (including infinitely near points) with multiplicity at least that of C minus one (see 
Section 4 for a more precise definition). One says that H € k[x, y] is an adjoint polynomial 
of F of degree n if it gives the dehomogenized equation of an adjoint curve of degree 
n. Adjoints may be computed by linear algebra from the resolution of singularities. One 
denotes by 

A C k[y] 

the vector subspace spanned by the remainders modulo (x) of adjoint polynomials of F 
of degree d — 2. Our main results assert that one can compute quickly both the rational 
and absolute factorizations of F from the knowledge of a basis of A. 

One assumes that fast Fourier transform can be used for polynomial multiplication, so 
that two univariate polynomials over k of degree < m can be multiplied in softly linear 
time 0(m). One denotes by 2 < uj < 3 the matrix multiplication complexity exponent. 

Theorem 1. There exists a deterministic algorithm that, given F satisfying (H) and 
given a basis of A, computes the rational factorization of F with one factorization in k[y] 
of degree d plus 

0(d 2 (d- s)"- 2 ) c 0{d u ) 
arithmetic operations over k, with s the number of irreducible rational factors of F. 

The vector space A contains enough information to compute too the absolute factoriza- 
tion of F. By absolute factorization, one means here the computation of a family of pairs 
of polynomials 

{(<3i,<7i), • • • , {Qr,q r )} 
where € k[t] is monic, Qi £ k[x,y,t] and where 

r 

i=l qi(a)=0 

is the irreducible decomposition of F over k. Note that such a representation is non 
unique. One obtains the following result : 

Theorem 2. Suppose that k has characteristic or greater than d(d — 1). There exists 
a deterministic algorithm that, given F satisfying hypothesis (H) and given a basis of A, 
computes the absolute factorization of F within 

0(d 2 {d - s)"- 2 + sd 3 ) C 0(d 4 ) 

arithmetic operations over k, with s the number of irreducible absolute factors. 
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Note that in contrast to the rational case, no univariate factorization is required for 
absolute factorization. Following (Cheze-Lecerf 2007), one obtains too a probabilistic 
approach in terms of computation trees. 



Theorem 3. Suppose that k has characteristic or greater than d(d — 1). Given F 
satisfying hypothesis (H) and given a basis of A, there exists a polynomial S G k[ti, . . . , 
of degree at most d(d— 1) and a family of computation trees parametrized by c £ k d such 
that 

• Any executable tree returns the absolute factorization of F ; 

• A tree is executable whenever S(c) ^ 0. 

The maximal cost of the trees is bounded by 

6{d 2 {d-sy~ 2 + d^) C O(d^) 

arithmetic operations over k. If the cardinality of k is infinite, the algorithm returns the 
correct answer with probability one. 



Up to our knowledge, all our complexities are smaller or equal to the actual best com- 
plexities for factorization of dense bivar iate polynomials that are known to be 
for rational deterministic (Lecerf 2007), 0(d 4 ) for absolute deterministic and 0(d 3 ) for 
absolute probabilistic (Cheze-Lecerf 20071. Thus, our results lead immediately to the 
following question : 

Question: Can we compute a basis of A fast enough for our method being useful for 
rational or absolute factorizations ? 



One obtains the following encouraging result, as a consequence of the Ricmann-Roch 
theorem for reducible curves. 



Theorem 4. Given a basis of the vector space of all adjoint polynomials of degree d— 2 , 
one computes a basis of A within 

0(d(g + d- s)(d - s)"- 2 ) C 0{d u+1 ) 

arithmetic operations over k, where g is the geometric genus of C. 



Note that there exist efficient algorithms for computing adjoint polynomials, by using 



Newton-Puiseux expansions ( Stadelmeyer- Winkler 1997) or integral basis (Mnuk 1997 



Deconinck-Van Hoeij 2001), but whose complexities have not been analyzed yet. Un- 



fortunately, it is a priori hopeless that Theorem 4 answers positively to our question : 
one needs the all resolution of singularities to compute adjoints and one expects factor- 
ization to be one step in desingularization (up to our knowledge, if k = ¥ p with p > d, 
the best complexity for computing all singular Puiseux expansions of F is 0(d 5 ) opera- 
tions ( Poteaux-Rybowicz 2011 )). However, one will show (Section 8) that it's enough to 
separate all local branches C in order to compute A. In particular, one needs not to desin- 
gularize irreducible branches. Moreover, one will see (Subsection 8) that in some cases, 



3 



our method adapts too to the case F(0, y) non separable, in which case one can use the 
combinatorial information given by the resolution of singularities along the line x = 
to speed-up the algorithm. In that spirit, the author recently developed a factorization 



algorithm based on the toric resolution of the singularities at infinity (Weimann 20101 



running in polynomial time in the volume of the Newton polytope, improving (Lecerf 
[2007] ) for sparse enough polynomials. 

To summarize, the method developed here may be useful for some special type of 
polynomials (smooth components, high singularities along a line, etc.) or more generally 
if one is given some extra input data concerning the singularities. Anyway, our results 
clarify the relations between normalization and factorization, giving a good point of 
view for comparing the classes of complexity of both operations and of various related 
algorithms (Newton-Puiseux, Hensel lifting, integral closure, etc.). It has to be noticed 
too that the strength of our approach depends strongly on further improvements in the 
algorithmic theory of singularities, especially on the Newton-Puiseux algorithm. 

The proofs of our main results rely on the structural sheaf sequence of a divisor on 
the normalization of C, combined with the Serre duality and with the residue theorem. 



Roughly speaking, our algorithms combine ideas developed in (Duval 1991 ) a nd (|Ragot| 
1997) (computing locally constant rational functions) with ideas in (Lecerf 2007) and 



( Cheze-Lecerf] 2007) (lifting and recombination of modular factors). Namely, one shows 
that one can recombine the factors modulo (x) from A without using Hensel lifting. In 
fact, one proves that absolute recombinations and adjoints modulo (x) determine each 
other by solving a d x d linear system over k (see Corollary 18 for a precise statement). 

Organization. One introduces the recombination problem and its relation to locally 
constant functions in Section 2. In Section 3, one proves our key result that gives con- 
ditions for lifting locally constant functions using residue theory and cohomology. In 
Section 4, one establishes the relation with adjoint polynomials and proves Theorem 4. 
One solves recombinations in Section 5 from which follow the proofs of Theorem 1, 2 and 
3 in Section 6. In Section 7, one discusses the case F(0,y) non separable and one illus- 
trates our approach on a simple example. In Section 8, one shows that the computation 
of A does not require the all resolution of singularities. Finally, one concludes in the last 
Section 9. 



2. Recombinations and locally constant functions. 



Our algorithms are related to ( Lecerf 2007 1 and ( Cheze-Lecerf 2007 ) , both methods 



being based on the recombination problem of the modulo (x) factors. One first explains 
this problem and then one relates it to the sheaves of locally constant functions on the 
normalizing curve. One keeps the same notations and hypothesis as in the introduction. 

2.1. Recombinations problems 

Let us consider the respective factorizations 

\f(x, y) = Fi(x, y) ■ ■ ■ F s (x, y) 
\F(0,y) = f 1 (y)---f n (y) 
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of F and F modulo (x) over k (recall that F(0,y) is assumed to be separable). Solving 
rational recombinations consists in computing the vectors 

^•) = (^\...,^))e{0,ir 

induced by the relations 

71 

F j (0,y) = l[f i (yyy\ j = l,...,s. 

i=l 

In the same way, let 

[F{x,y) = F 1 {x,y)---F- s {x,y) 
\F(0,y) = f 1 (y)---f d (y) 

be the respective factorizations of F and F modulo (x) over k. Solving absolute recom- 
binations consists in computing the vectors 

u^ = (p{ j \...,p^)e{o,i} d 

induced by the relations 

d 

F,(P,V) = t[Mv) D ' \ 3 = I,..., s. 

i=i 

The following picture illustrates the absolute recombinations when C is union of a cubic 
and a conic. 



i/W = (1,0,1,0,1), i>< 2 >= (0,1,0,1,0) 

In this article we mainly pay attention to the recombination problems, the irreducible 
factorization of F then following with a fast multi-factor Hensel lifting (combined with 
a partial fraction decomposition algorithm in the absolute case). The main idea is to 
interpret the recombination problem as a cohomological problem of lifting sections. 



2.2. Solving recombinations via lifting sections 

All schemes and properties (connectivity, irreducibility) are considered over the base 
field k. One may think a point over a scheme X over a A: as a collection of points in the 
extension X — X <g)fc k that are conjugated under the Galois group of k/k. 

Let C and C be the respective Zariski closures of the affine curves F = and x — to 
the projective plane P 2 . Let 

7T : X -> P 2 

be the standard embedded resolution of C. One denotes by C and L the respective strict 
transforms of C and C by tt. The inclusion of the zero-dimensional subscheme 

Z :=CnL 
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into C induces a restriction morphism 

a : H°(O c ) ^ H°(O z ) 

between the respective fc-vector spaces of regular functions on C and Z. Note that both 
vector spaces may be identified with the sets of locally constant functions on C and Z. 
The map a is injective since Z has at least one point on each component of C. The 
following two subsections are dedicated to show that the computation of the cokernel of 
a permits to solve both the rational and the absolute recombination problems. 



2.2.1. The rational case. 

Since F(0, y) has degree d, Z is an affine zero-dimensional subscheme whose ring of 
regular functions may be identified with the finite fc-algebra 

H {Gz) ~ m = Jmwy (1) 

Since F(Q, y) is separable, its rational factorization induces an isomorphism 

Thus Z has n connected components (closed points) pi,...,p n corresponding to the 
maximal ideals of the ring H°(Oz) generated by the fi's. The natural inclusions 

k\y] 

combined with (1) and (2) induce the inclusion 

k n being identified with the subspace of locally constant functions on Z that take value in 
k, that is (vi, . . . , v n ) e k n sends p, to Vi (in general, a function on Z takes values in the 
various residue fields k[y]/(fi)). The map a introduced before is related to recombinations 
by the following lemma: 

Lemma 5. The vector subspace W C k n defined by 

W := k n n Im(a) 
admits (u^\ . . . , ^^ s ') as reduced echelon basis (up to reordering). 



Proof. By definition, v e W if and only if it's the restriction to Z of a locally constant 
fc-valued function on C. Since C is smooth, it has s connected components C\ . . . , C s 
corresponding to the prime rational factors of F . Thus v G W if and only if v is k- 
valued and constant along Cj n L for j — 1, . . . , s. One deduces that dim^ W — s and 
that e W for j = 1, . . . , s. Since the i/^'s have {0, l}-coordinates and are pairwise 
orthogonal vectors in k n , they form up to reordering the reduced echelon basis of W. □ 

By Lemma 5, the recombination problem over k is reduced to compute first the rational 
factorization of F(0,y) (inducing the inclusion k n C H°(Oz)), and then the cokernel of 
a. 
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2.2.2. The absolute case. 



The relations between locally constant functions and absolute factorization is explored 
in (Duval 1991 ) where the author determines one absolute factor a time from a basis of 
the regular functions on C Xj I. One rather relates here regular functions on C to the 
recombination algorithm in ( Cheze-Lecerf 2007 ) and one computes all irreducible factors 
simultaneously by using multi-factor Hensel lifting. One first proves : 



Lemma 6. One has equalities dirrikH (Oz) = d and dim k H°(Oc) = s. 



Proof. First equality is clear from (1). Since H a (Oc) is a finite dimensional fc-vector 
space, one has 

dim k H°(O c ) = dim- k H°(O c ) ® k k 
Let C :— C Xj. k be the geometrical scheme ass ociated to C by extending the base field 



k to its algebraic closure k. One has (Liu, 2002 prop. 1.24 p. 85) 



H\O c ) ® k k = H°(O e ). 

Since C is smooth, it's the disjoint union of s irreducible components C\ . . . , C s corre- 
sponding in an obvious way to the prime absolute factors of F. It follows that one has 
an isomorphism of k- vector spaces 

Since H°(Oq.) = k, one has dim k H°(Oc) = s so that dim k H°(O c ) = s. □ 



Let fa,..., <pd be the roots of F(0,y) in k. The identification (1) gives rise to the 
multi-evaluation isomorphism 

ev: H°{O z )®k~k ^> k d 

v i — > v := {v(<pi), . . .,v{<j>d))- 

The next lemma shows that solving absolute recombinations reduces to compute Im(a) 
and to apply the evaluation map ev. One endows k d with its canonical basis. 

Lemma 7. The vector subspace W C k d defined by 

W := ev(Im(a) ®fc k) 
admits (f'*- 1 -', . . . , i>^) as reduced echelon basis (up to reordering). 



Proof. Let Z := Z x k k. The map ev induces an identification 

k d = H°(O z ), 

where v = (pi, . . . ,i>d) G k d is identified with the locally constant function that sends 
each closed point pi € Z to vi. Since Z contains at least one point of each connected 
component of C, the restriction map 

a : H°(O c ) ^ H°(0 2 ) 



7 



is injective. By definition, W = Im(a) so that dim^W = dim^H \0 'g) = s by the 
proof of Lemma 6. Each vector i/Ci) being constant on C\ n L, . . . , C s (1 L, it extends to 
a function on (7. So € for j = 1, . . . , s. Since the i/^'s have {0, l}-coordinates 
and are pairwise orthogonal in they form up to reordering the reduced echelon basis 
o£W. □ 



3. Lifting sections using residues 



The previous section shows that recombinations may be reduced to compute the cok- 
ernel of the restriction morphism 

a : H°(O c ) ^ H°(O z ). 



To this aim, one introduces residues. There is an extensive litterature ( Couvreur 2009 



Serre 1988 Lipman 2011 Tate 1968 Vakil 20081 concerning residues on curves and 



surfaces. 

Let ujc be the sheaf of regular differential 1-forms over C (the dualizing sheaf) and let 
ljc{Z) be the sheaf of meromorphic 1-forms with polar divisor bounded by Z. Let jjgC 
with residue field k p and let ip £ u>c. P {Z) be a germ of meromorphic form at p. For any 
uniformizcr t of C at p, there exists a unique formal series h G such that 



h(t)dt 



One defines the residue of ip at p as 



resp ip := Tr p [ft(0)] 5 

where Tr p : k p — > fc is the trace map. This definition does not depend on the choice of 



the uniformizer (see for instance (Serre 1988)). The map res p is k- linear and vanishes 
on regular forms. In particular, if v £ Oz, P has a local lifting v to Oc, P , one checks that 
the definition 

res p (i/ , 0) := Tes p (vip) 
does not depend on the choice of the lifting. One obtains the following key result. 



Proposition 8. There is an exact sequence of k-vector spaces 

— > H°(O c ) H°(O z ) A H°(ujc(Z)y -A H°{ujcY — > 
where v stands for the dual and where R associates to v the linear form 



res Pi (utp). 

i=l 

In particular, dim H°(ujc{Z)) = g + d — s where g is the geometric genus of C (sum of 
the genus of the irreducible components). 
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Proof. Let ujz be the dualizing sheaf of Z. One has loz — Hom(Oz,k), and the local 
duality theorem gives a short exact sequence (the adjunction formula) 

— > u>c — > Lo c (Z) ^ u>z — > 0, 
where the residue map Res is locally defined on an open set U C C as 

ResuW : O z (U) — > k 

v i — > res p (vip). 

The associated long exact cohomology sequence is 

-+ ff°(wc) fl°(w c (Z)) ^ s H°(w z ) -+ H\w c ) H\u c {Z)). (4) 

By the duality of Serre, one has isomorphisms 

H\uj c )^H {O c ) y and H\uj c {Z)) ^ H°{O c {-Z)) y =Q, 

the last vanishing property because Z as at least one point on each connected component 
of C. The dual sequence of (4) becomes 

-> ff°(O c ) 4 #°(0z) 4 H°{io c {Z)) y 4 H°(lj c ) V -> 

where i? is dual to i?es, that is 

This shows the exact sequence of Proposition 8. This sequence induces equality 
h°(u c (Z)) = h {uj c ) + h°(O z ) - h°(O c ) = g + d~s, 

last equality using Lemma 6 and using that h°(wc) coincides with the geometric genus 
of C. This ends the proof. □ 



Remark 9. The inclusion Im(a) C ker(R) follows from the residue theorem that asserts 
that 

res p ip = 

for all connected component Cj of C and all rational 1-form ip on C. 



Remark 10. In the case C irreducible over k, the equality 

dimH°{u) C {Z)) = g + d-l 
given by Proposition 8 follows from the theorem of Riemann-Roch for curves. 
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4. Relations with adjoint polynomials 



One relates now holomorphic forms with adjoint polynomials. One denotes by S the 
set of singular points of C, including all infinitely near points. For each p£S, there is a 
decomposition of ir 

x^hx p ^x p ^v 2 

such that p is a closed point of the intermediary surface X p and tt p is the blow-up at p. 
Let E p be the exceptional divisor of tt p and E p its total transform under tti . One denotes 
by m p the multiplicity at p of the strict transform of C under the map 7T2 . 

Definition 11. The adjoint divisor of F is the exceptional effective divisor 

E:=^2(mp-1)E P . 

pes 

An adjoint curve of C is an effective divisor J)cP 2 that satisfies 

7T*(L>) > E. 

An adjoint polynomial of F of degree < m is a polynomial giving the dehomogeneised 
affine equation of an adjoint curve of degree to. 

In other words, adjoints of F are those polynomials vanishing at the singular points 
of C with high enough multiplicities. Adjoints carry out precious informations about the 
geometry of C. In particular, it is well known that they are deeply related to the sheaf 
ujc of regular forms on the normalized curve. Let us denote by 

Adj(m) C k[x, y] 

the fc-vector subspace generated by adjoint polynomials of F of degree < to. One has the 
following proposition: 

Proposition 12. For all integers m < 2, one has an isomorphism 



Adj(d - 3 + to) H°(u c (mZ)) 
/ Hdx \ 
n \x m d„F) 



H 



\c 



Proof. In order to relate adjoints with differential forms, one introduces the conductor 

Ac := Hom 0c (ir*Oc,O c ) 

of the normalization of C. It is an ideal sheaf of Oc, related to the dualizing sheaf u>c of 
C by the formula 



(see for instance Szpiro| 1979 p. 25). Since the morphism tt : C — >• C is affine one has 



H°(C,u c (mZ)) =H°(C,n.(LJ C (mZ))) = H° {C ,Lo c {mC) ® Ac) , (5) 
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last equality following from the projection formula (recall that Z = COL and ir*£ = L). 
Let A be the inverse ideal sheaf of Ac under the restriction Op2 —> Oc- One has the 
short exact sequence 







A — > Ar 



0. 



(6) 



Tensoring (6) with the invertible sheaf ttp 2 {C + raL), and using the adjunction formula, 
one obtains the exact sequence 







0. 



fip 2 (m£) — > njjL (C + mC) ®A cu c (m£) <8> A c 

Here, RP is the Poincare residue map, defined outside the singular locus of C as 

' Hdx Ady\ / Hdx \ 
Fx m ) ~ \x m d v FJ\c' 



RP 



(7) 
(8) 



For m < 2 and i = 0, 1, one has ( |Hartshorne] |1977| Theorem 5.1 p. 225) 



,fi^(m£)) =if i (P 2 ,0 P 2(m-3)) = O 



which, combined with the long exact cohomological sequence of (7), gives an isomorphism 
RP : H°{¥ 2 ,ttl 2 {C + mC)®A) H a {C , LJ C {mC) ® Ac) ■ (9) 
By (Szpiro 1979 Proposition p. 33), one has an isomorphisrrp"] 

Adj{d - 3 + m) ^ H°{P 2 ,nl 2 (C + mC) ® A) 

Hdr. , 

H 



Hdx A dy 



Fx m 

which combined with (5), (8) and (9) gives the isomorphism of Proposition 12. □ 



Remark 13. The isomorphism of Proposition 12 for m = and C irreducible is known 



as the Gorenstein Theorem (Gorenstein 1952) which asserts that the adjoint curves of 



degree d — 3 cut out on C the complete canonical system of its normalization. For a nice 



down-to-earth presentation of adjoints and conductors, see (Fulton 2004 ) . 



Recall from the introduction that one defines Ac k[y] to be the image of the projection 

Adj(d-2)^k[y] 

H^H(0,y). 

Corollary 14. One has equality dira^ A = d ~ s. 



One checks that the irreducibility assumption made in ( Szpiro 
since the proof is local. 



1979 



Proposition p. 33) can be removed 
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Proof. If H £ Adj(d — 2) satisfies H(0,y) = 0, then H(x,y) = xH'(x,y) for some poly- 
nomial H' . Since the line x = does not contain any singularities of C, H' is necessarily 
an adjoint of F of degree d — 3. In other words, one has an exact sequence of fc-vector 
spaces 

— > Adj(d - 3) — > Adj(d - 2) — > A — > (10) 

where the first map is the injective "multiplication by x n map and the second map is the 
restriction to x = 0. It follows that 

dim(A) = dim Adj(d - 2) - dim Adj(d - 3) 

= h°(u c (Z))-h (uc) 
= d— s, 

second equality using Proposition 12 and last equality using Proposition 8. □ 
The proof of Theorem 4 follows. 

Corollary 15. (Proof of Theorem 3). Given a basis of Adj(d— 2), one can compute a 
basis of A within 

0({d - l){g + d - s)(d - sT" 2 ) C 0((T +1 ) 
arithmetic operations over k. 



Proof. Consider the matrix N whose set of rows is a basis of Adj(d — 2) evaluated at 
x = 0, expressed in the natural basis of k[y]. So N has d — 1 columns and g + d — s 
rows (use Propositions 8 and 12). By (10), one has A — Im(N), and a basis of A can be 
computed within the expected complexity ( Storjohann| |2000| Theorem 2.10). The upper 
bound 0(d u+1 ) follows from the well known inequality g < (d — l)(d — 2)/2. □ 



5. Recombinations follow 

One has now all necessary information for solving recombinations. Let us consider first 
the rational case. 

Corollary 16. One has an exact sequence of k -vector spaces 

0^ (v (1 \...,v (s ^) — >k n A A v 
where T sends v — [y\, . . . , v n ) to the linear map 

1=1 \/i(0)=o y ' 1 
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Proof. By Lemma 5, one has 

=Im(a)nk n 

where one identifies k n C H°(Oz) with the subspace of locally constant k- valued func- 
tions on Z. Proposition 8 induces equality 

n 

Im(a) n k n = {1/ e fc n , 5Zi/ i rea P4 (^) = V V G fl°(w c (^)}. 

i=i 

Let us compute the involved residues. By Proposition 12, ip e H°(ljc{Z)) is equal to 

for a unique € Adj(d — 2). Let 0c,pj be the completion of the regular local ring 0c,p; 
with respect to its maximal ideal associated to p^. The residue field of C at Pi is equal to 

k -*M 
Pi " (/,)" 

The map 7r being an isomorphism in a neighborhood of Pi, one has an isomorphism 

dc, Pt ^k Pi [[t]\ 

7T X I > t 

n*y i — > a(t) 

where a € fc Pi [[i]] is the unique series such that a(0) is the residue class of y in k p . and 
F(t, a(t)) = 0. In such a local system of coordinates, ip is equal to 

_ H(t,a(t)) dt 
d y F(t,a(t))~t 
and it follows from the definition of residues that 

res (*) Tr ( J^lA) ST ( JUML) 
" W ~ " V d y F(0, a(0))) f ^ =Q \ d y F(0, $ ) ■ 

Corollary 16 follows. □ 

Remark 17. One always has (1, . . . , 1) G ker(T). This is nothing else than the Lagrange 
interpolation formula. 

Let us now consider the absolute case. 

Corollary 18. One has an exact sequence ofk-vector spaces 

0^ (i> (1) ,...,P (5) ) — ^k d ^ A v (g> k k^0 
where T sends v — (D\, . . . , D4) to the linear form 

H{4>i) 
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Proof. Apply Proposition 8 and repeat the proof of Corollary 16 over k, with the curve 
C replacing C. Surjectivity of T follows from Corollary 14. □ 



Remark 19. In ( Cheze-Lecerf 2007 Lecerf 20071, the authors solve recombinations 



using a system of 0(d 2 ) equations. Corollary 16 and Corollary 18 give a much smaller 
number d— s of equations for recombinations. Moreover, the map T being surjective, d—s 
is the expected minimal number of linear conditions for recombinations in the absolute 



6. Proofs of Theorems 1, 2 and 3. 

In all what follows, one assumes that fast Fourier transform is used for polynomial 
multiplication, so that two univariate polynomials over k of degree < m can be multiplied 
in softly linear time 0(m). 

6.1. Proof of Theorem 1. 

One obtains the following algorithm. 

Algorithm 1 (deterministic rational factorization) 

Input : F € k[x, y] that satisfies hypothesis (H). 
Output : The rational factorization of F. 

• Step 1. Compute a basis of A. 

• Step 2. If dim A = 1, F is irreducible. Otherwise, compute the irreducible factors 
fl, ■■-,fn of F(0,y) over k. 

• Step 3. If n = 1, F is irreducible. Otherwise, build the matrix M of the map T of 
Corollary 16 by using Newton identities. 

• Step 4. Compute the reduced echelon normal basis of ker(M). One obtains the recom- 
bination vectors v^ x \ . . . , v^ s \ 

• Step 5. Compute the factorization of -F(0, y) induced by the recombination vectors and 
lift it to the rational factorization of F. 



Proposition 20. (Proof of Theorem 1.) Algorithm 1 is deterministic and correct. Steps 
3, 4 and 5 take at most 

0(n(d - s)(d - s)"" 2 + d 2 ) C 0{d u ) 
arithmetic operations over k. 



Proof. The algorithm is deterministic and correct thanks to Corollary 16. Let us describe 
in more details the content and the complexity of steps 3 to 5. 
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Step 3. In order to build the matrix M, one has to compute 

H(y) 



Tr. 



: 'Ad v F(0,y)J 



Tn 



ajTr k , 



XV s )- 



(ii) 



for all i = 1, . . . , n and for all H running a basis of A. Inversion of d y F(0, y) and multi- 
plication by H in k[y]/(fi) take 0{ni) operations in k. Then H/d y F(0,y) 6 k[y]/(fi) is 
uniquely represented as a polynomial a(y) = do + ■ ■ ■ + ci ni ^iy ni ~ l with coefficients in k 
and 

g(y) \ = r v 

Thanks to the Newton identities, one can compute recursively the trace of y° from the 
traces of smaller powers of y and from the coefficients of fi with j multiplications and 
j additions. So one computes traces of all involved powers of y within 0(n 2 ) operations 
over k. Given these traces, and using (11), one computes the trace of H/d y F with 2n^ 
operations for each H £ A. By Corollary 14, it follows that step 3 costs C( n i + 

2rii(d — s)) C 0(d 2 ) operations over k. 

Step 4. The matrix M has size (d — s) x n and rank d — s. One can compute the 
reduced echelon normal basis of the kernel of M within (D(n(d— s)(d— s)"~ 2 ) operations 
(( |Storjohann|[2000| , Theorem 2.10). 

Step 5. Given a vector — (i/,- ) G {0, 1}" of the reduced echelon basis, one com- 

u) 

putes Fj(0,y) — YlfXv)^ l° r each rational irreducible factor Fj of F. This requires 



O(deg(Fj(0,y))) operations by the sub-product tree technique (Lecerf 2007 proof of 
Prop. 6), so a total cost of O(d) operations. To compute the Fj's, it's now enough to 
lift the induced equality F(0, y) — i*i(0, y) ■ ■ ■ F s (0, y) modulo (x) up to precis ion mod 
ulo (x d+1 ). This costs 0(d 2 ) operations by using Newton quadratic iteration (Gathen- 
Gerhard[ [2003} Theorem 15.18). □ 



6.2. Proofs of Theorems 2 and 3 

In the absolute case, the delicate point is that Corollary 18 does not permit to solve 
recombinations with linear algebra over k. Moreover, it neither permits to describe the 
smallest finite extensions over which the irreducible absolute factors of F are defined. 
To solve this problem, one rather relies our approach with the algorithms 8 and 9 in 



( Cheze-Lecerf 2007 ) , where the authors use the absolute partial fraction decomposition 
algorithm of Lazard-Rioboo-Trager ( Lazard-Rioboo 1990[ ) . 

Let <f> be the residue class of y in the ring A := k[y]/(F(0, y)). Any element b £ A can 
be uniquely represented as a finite sum 

d-l 

where coeff(6, (j) 1 ) :— 6, belongs to k. One introduces 



a 

L:=$yV£k d , Vj coeff ( 



i=l 



g(0) 

d y F(o,<i>y 



VH£A\ 
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The vector space L is related to the absolute recombinations by the following lemma. 



Lemma 21. Let V be the Vandermonde matrix of the roots 4>i, . . . , 4>d of F(0, y). One 
has an isomorphism 

V 1 : (£«,..., z^)) -^L® k k. 
In particular, one has an isomorphism of k-vector spaces 

B : Im(a) L 

where B = (Trcj)' l ' +J )i i j=o,...,d-i, with Tr : A — > k the usual trace map. 



Proof. One follows the proof of Proposition 4 in (Cheze-Lecerf 2007). Let (v%, . . . , Vd) = 
V t (wi, . . . , Wd) and let b e A. One has 

d d d 

Vi coeff (b, = X] (Zl w 3$) coeff ( fe ' 

i—1 i—1 j — 1 

d d d 

= E w i (E coeff ( 6 ' = E 

j=i i=i i=i 

The first point then follows from Corollary 18 by taking b — H(4>) /d y F(0, 4>). The second 
point follows from Lemma 7 since V is the matrix of the evaluation map and B = VV 
is the matrix of traces. □ 

One can now rely on the factorization algorithms developed by Cheze-Lecerf in the abso- 



lute case. One refers to their article (Cheze-Lecerf, 2007) for details on the relations be- 



tween absolute recombinations, absolute partial fraction decomposition, absolute Hensel 
lifting and absolute factorization. 



Algorithm 2 (deterministic absolute factorization). 

Input : F e k[x,y] that satisfies hypothesis (H), with k a field of characteristic or 
greater than d(d — 1). 

Output : The absolute factorization of F. 

• Step 1. Compute a basis of A. 

• Step 2. Compute a basis of L. 



• Step 3. Call Algorithm 8 in (Cheze-Lecerf 20071 with input F and the basis of L. 



Proposition 22. (Proof of Theorem 2.) Algorithm 2 is deterministic and correct. Steps 
2 and 3 take at most 

didid-s)"- 1 + sd 3 ) cd(d 4 ) 

arithmetic operations over k. 



1G 



Proof. The algorithm is correct thanks to Lemma 21 combined with Proposition 4 p. 15 



and Theorem 5 p. 15 in (Cheze-Lecerf 2007). By definition, one has L = ker(N), where 
the matrix N is built from a basis of A using one inversion in A and (d—s) multiplications 
in A, so 0(d(d— s)) operations over k. Then, computing a basis of L — ker(N) requires 
(D(d(d— s) w_1 ) operations over k. Finally, step 3 costs 0(sd 3 ) operations over k thanks 
to Proposition 10 p. 24 in ( |Cheze-Lecerf] [2007| . □ 



The cost of Algorithm 2 is dominated by the separation of residues in Algorithm 8 of 



(Cheze-Lecerf 2007) that ensures that the call to the Lazard-Rioboo-Trager algorithm 



returns a correct answer. If one rather deals with a random linear combination of the 
vectors of a basis of L, one obtains a probabilistic algorithm with smaller complexity. 



Algorithm 3 (probabilistic absolute factorization). 

Input : F E k[x,y] that satisfies hypothesis (H), with k a field of characteristic or 
greater than d(d — 1). 

Output : The absolute factorization of F. 

• Step 1. Compute a basis of A. 

• Step 2. Compute a basis of L. 

• Step 3. Choose c £ k s , where s = d — dim(A). 



Step 4. Call Algorithm 9 in (Cheze-Lecerf 20071 with input F, the basis of L and c. 



Proposition 23. (Proof of Theorem 3.) Algorithm 3 either stops prematurely or return 
a correct answer. Moreover, there exists a polynomial S € k[C\, . . . , C s ] of degree at most 
s(s — 1) such that the answer is correct whenever S(c) ^ 0. In any cases, steps 2, 3 and 
4 take at most 



0(d(d-s) a 



arithmetic operations over k. 



Proof. The proof is the same as for Proposition 22, using now Proposition 11 p. 25 in 
(ICheze-Lecerfl [20071). □ 



Remark 24. The restriction hypothesis on the characteristic of k ensures the possibility 



to separate the residues and to apply a fast absolute multi- factor Hensel lifting in ( Cheze- 



Lecerf, 2007 Algorithm 9). 



7. The case .F(0, y) non separable 

When F(0,y) is not separable modulo (x), one is tempted to choose another fiber 
x = a for which F(a,y) that satisfies hypothesis (H). There are two main reasons to 
develop a recombination algorithm along a critical fiber. First, when the field k has small 



17 



positive characteristic, a regular fiber may not exist. Second, one shows here that working 
along a singular fiber may in fact be an opportunity to speed-up the algorithm. 

In order to simplify, one supposes here that k — k. One supposes too that F(Q,y) has 
degree d (the general case follows easily by computing residues at y = oo). 

Our results generalise well to the non separable case, the main difference being related 
to the computation of residues. Let Z — C fl L where L — tt*(C). In contrast to the 
previous sections, L and Z need not to be reduced anymore. The support of Z consists 
now in r < d closed points pi , . . . , p r in one-to-one correspondance with the irreducible 
analytic branches of C along the line £. The recombination vectors may now be defined 
in the smaller ambient space k r , where fjP' G k r is defined to have ith coordinate equal 
to 1 if pi <G Cj and equal to otherwise. By identifying k r with the vector subspace of 
H°(Oz) of locally constant functions on Z with values in k and with zero nilpotent part, 
one obtains the analoguous of Lemma 5 

(i/W,...,^) =Im(a)nk r , (12) 

where a : H°(Oc) — > H°(Oz) still stands for the restriction map. One obtains the 
following generalization of Corollary 16 : 

Proposition 25. One has an exact sequence of k-vector spaces 

—►<!/«,..., i/W) -^k r ^>A V 
where T sends v — [y\, . . . , v r ) to the linear map 



Proof. The exact sequence of Proposition 8 remains valid in this new context. Combined 
with (12), one obtains that (z/ 1 ), . . . , z/ s )) = kerf where f : k r -> H°{lo c {Z)Y sends 
v = (vi, . . . , v r ) to the linear map 



The main difference concerns the computation of residues. Let if> € H (wc(Z)). By the 
proof of Proposition 12, one has 



i\> = RP C (V), 



* = 7T* 



H (x, y)dx A dy 
xF(x,y) 



e H°(n 2 x (c + L)) 



for a unique H 6 Adj(d — 2), and where RPc stands for the Poincare residue along C . 
Let p 6 C fl L. By ( Couvreur] 2009, Theorem 6.3 and Remark 6.9), one obtains equality 



res p (4>) = res p (RPc(V))) = res p (i?P L (*)) 
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where the last residue stands for residue of 1-form on L. Since L = it* (£) and the Poincare 
residue commutes with the pull-back, one obtains equality 

DD ( J ' H(x,y)dxAdy \\ J ( H(x,y)dx A dy \\ J H(0,y)dy \ 
RPL { n { xF{ X ,y) )) =7T [ RPC { xF(x,y) ) ) = * J " 

Proposition 25 follows. □ 

So as soon as the curve C has a small number of irreducible branches intersecting C, 
Proposition 25 permits to solve the recombination problem in a smaller ambient space. 
The price to pay is that one can a priori not compute residues directly in P 2 as in Section 
2, but one may really need to compute residues in X, using local coordinates or Puiseux 
series. Nevertheless, in the important case of C locally irreducible at (0, y p ) e Cn£, there 
is exactly one point p G C such that ir(p) — (0,y p ) and the residue can be computed 
directly in P 2 : 

r€Sp y \F%y))) reSy A F (0,y)) ■ 

Of course y p may be now be a multiple root of F(0,y) so the residue computation may 
involve higher order derivatives of H and F. 

Example. Let us illustrate Proposition 25 on a simple example. Suppose that one wants 
to factorize 

F(x, y) = y 5 + y 4 - xy 3 - y 3 - 2xy 2 - y 2 + x 2 + xy + x. 

overQ. The irreducible factorization of F mod (x) is 

F(Q,y)=y 2 (y + l) 2 (y-l), 

with two double roots —1 and and one simple root 1. Since d x F(0, —1) and d x F(0, 0) do 
not vanish, the curve C is smooth and tangent to C at these two points, and transversal 
to C at (0, 1). In particular, C has only 3 irreducible branches intersecting C at distinct 
points and the recombinations will hold in the ambient space k 3 rather than in the bigger 
space k 5 inherent to a choice of a regular fiber. 

One has here 

Adj(d - 2) = (y 3 - y - 1, y 2 - x, y 2 - x, xy 2 - x 2 , y 3 - xy) 
from which it follows that 

A=(y + l,y 2 ,y 3 ). 

By Lemma 6, the curve C has deg(F) — dim(A) = 2 absolute irreducible components. Let 
H e A. Since F is locally irreducible at (0,0), one has equality 

( H{y)dy \ ( H(y)dy \ 

reS °l F(0^j) - reS ° W(y + iny-l) ) =H{ ° ) + H{0) - 

In the same way, a simple calculation gives 

( ^ = -2^(-i)- 5g( -i) and nai (^m = m. 

l \F(0,y)J 4 yF(0,y)J 4 
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One deduces that the 3x3 matrix of the map T in Proposition 25 is 

f -l/2 -1/4 -1/4^ 
M = 

v 1/2 1/4 1/4 j 

so that ker(M) = ((0, 1, 0), (1, 0, 1)). One deduces the irreducible rational decomposition 

C = Ci U C 2 

where C\ is a conic tangent to C at (0,0) and C 2 is a, (smooth) cubic tangent to C at 
(0, 1) and transversal to C at (0, —1). 



x=0 



Since the induced factors of F are coprime modulo (x), one can compute them use a 
multifactor Hensel lifting (Gathen- Gerhard] \2005\ Algorithm 15.17) up to a sufficiently 
high precision (mod (x 3 ) in our case). One obtains finally 

F(x,y) = (y 2 -x)((y + l) 2 (y-l)-x). 



In summarize, one has shown here that working along a critical fiber may be an 
opportunity to speed-up the algorithm, at least when F satisfies the weaker hypothesis 

(H') C is analytically irreducible at each point ol C n C. 

First, the univariate factorization of F(0, y) is faster since it is reduced to a fast separable 
factorization (Lecerf 20081 plus some univariate factorizations of smaller degrees. Second, 



recombinations are faster since they hold in a smaller ambient space of dimension the 
number of distinct roots of F(0,y) (or irreducible factors in the rational case). This fact 
is well illustrated in a previous work of the author (Weimann 20101 who developed a 



lifting and recombination algorithm based on the toric resolution of the singularities of 
C along the line at infinity. 



Remark 26. One can show that under hypothesis (H'), building the matrix of the 
map T has the same cost in the separable and non separable cases. The fact that the 
computations of the residues may involve /i-adic expansions with higher precision for each 
irreducible factors fa of F(0,y) is compensed by the fact that the sum of the degrees of 
the fa's decreases. 



Remark 27. Although Proposition 25 still permit to solve recombinations even when 
F does not satisfy (H'), the problem resides in the fact that the irreducible factors of 
F may not be coprime modulo (x) and can not be computed with Hensel's lemma. This 
problem will be explored in a further work. 
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8. Don't touch the cusps 

It turns out that the computation of A does not necessarily require to compute the 
all resolution of singularities of C. Namely, let us consider the factorization of it 

X — > X Q ^ P 2 

where ttq is the minimal composition of blow-ups under which the strict transform Co of 
C is every where locally irreducible. Then one can check that all our results (Lemmas 5, 6, 
7 and the key Proposition 12) remain valid with Co replacing C and with the arithmetic 
genus p a (Co) > g of Co replacing the geometric genus of C (the proofs mainly only use 
that the irreducible and connected components of C coincide). Then, one checks easily 
that one has an exact sequence 

— > Adj (d - 3) ^4 Adj (d - 2) ^5 A — > 

where Adjo(k) is defined similarly as Adj(k) with the map itq replacing n in Definition 
11. For instance, one needs not to desmgularize cusps for computing A, which is of course 
natural from our factorization point of view. Note that there are easy local irreducibil- 
ity sufficient criterions that can be directly read off from the Newton polygon of the 
singularity (for cusps for instance). 



9. Conclusion 



One has established the bridge between locally constant functions (Duval 1991 Ragot 



1997) and lifting and recombinations algorithms ( |Lecerf 2007 Cheze-Lecerf 2007). One 
has shown that the computation of adjoint polynomials allows to solve the recombinations 
problems without using Hensel lifting and with the expected number of linear equations. 
Although one believes that our approach uses too strongly the geometry of singularities, it 
may be useful for some particular polynomials, for instance if the irreducible components 
are smooth and intersect in few points. Moreover, one has discussed the possibility of 
speed-up the algorithm when F is not separable modulo (x). In a further work, one 
plans to develop some intermediary algorithms based only on the resolution of some 



of the singularities, in the vein of (Weimann 2010) that uses the toric resolution of 



the singularities at infinity. In any cases, the power of using singularities for factorization 
depends strongly on complexity issues in the algorithmic theory of singularities, especially 
on the Newton-Puiseux algorithm. 
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